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ABSTRACT 


In the past decade, electronic devices based on graphene and the related two-dimensional 
(2D) materials have exhibited outstanding figures of merit. However, so far, the fabrication of 
two-dimensional diodes, as elementary building blocks of electronic devices, still relies on 
manual or semi-automated handling processing. To unleash their commercial potential, the 
integration of 2D materials into a fully-automated fabrication line is a critical step. Here, we 
elucidate the focused ion-beam writing as an automated approach to construct lateral diodes 
on a 2D heterostructure (MoSe,/G) consisting of the stacked monolayer graphene and MoSe;. 
Se-defects generated by focused ion writing endow the 2D heterostructure with unique 
electronic properties like the adjustable work function and the quasi-metallic state, which 
allows for the construction of the barrier at the boundary of the writing and non-writing 
region. Benefiting from this feature, the ion-beam-written heterostructure is used to realize 
rectifying and current regulating diodes. Exhibiting comparable performance to traditional 
diodes, the rectifying diode has a rectification ratio of ~10*, while the current regulative diode 
has a dynamic resistance larger than 4.5 MQ. Furthermore, to illustrate practical applications 
of these diodes in digital logic electronics, AND and OR logic gates are directly inscribed on 
the heterostructure by ion beams. Our work demonstrates the focused ion-beam writing as an 


additional strategy for direct-writing of 2D diodes on graphene-based heterostructures. 


INTRODUCTION 


The diode is the rudimentary building block for electronic circuits [1] for applications of 
logic and filter circuits. With exceptional electronic properties, two-dimensional (2D) 
materials, such as graphene [2, 3] and transition metal dichalcogenides (TMDCs) [4], have 
skyrocketed as candidate materials for the future evolution of electronics, attracting immense 
interests among the electron-device community [5]. Towards this purpose, many approaches 
for the 2D diode fabrication have been intensively proposed, such as stacking 2D layers with 
dissimilar charge-carrier types [6-8], hetero-atom doping [9, 10], electrostatic doping, [11, 12] 
and chemical doping [13, 14]. Nevertheless, these approaches require complete or partial 
manual operations, considerably limiting their commercial applications. For example, the 
integration of graphene with semiconductors forms a graphene/semiconductor (G/S) Schottky 
junction, which is recognized as a promising candidate for real-life applications [15, 16]. 
Nevertheless, the multiple site-directed transfers of graphene onto the substrate, as the crucial 
step in the preparation of G/S junctions, lack the benefits of fully-automated manufacturing 
tools up to now [15, 16]. The manual or semi-automated handling of G/S junctions brings 
challenges for its widespread applications, such as device reproducibility, fabrication yield, 
and reliability [16]. Therefore, an approach with a fully-automated or well-established 
process-line would represent a significant leap forward for the applications of 2D electronic 
devices [15, 16]. 

Focused ion beam (FIB) is a highly-automated technique to micromachine the target 
surface with a high spatial resolution of ~10 nm [17, 18]. Commercial FIB systems are 
characterized by a Ga' ion beam with an energy of 30 keV. While being commonly used for 
milling at a high fluence of Ga” ions, the FIB process can also be employed at low fluence to 
give rise to several effects in 2D materials, such as the introduction of defects in selected 
regions [19]. Firstly, through controlling specific FIB parameters, the sputtering of target 


atoms by FIB processing can be limited to the top atomic layer; secondly, atoms are 
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specifically removed from the target, enabling the defect-pattern on the target surface; and 
thirdly, defect-engineering (DE) triggered by FIB bombardment is capable of tailoring the 
electronic properties of materials for desired applications. The above effects are incredibly 
crucial and attractive for 2D transition metal dichalcogenides (e.g., MoSe;), which contain 
three layers of atoms. By using FIB, defects (1.e., Se-vacancies through sputtered Se atoms) 
can be introduced locally in the first atomic layer to achieve modification of the electronic 
structure in a specific micro-region [19-21]. Furthermore, through stacking the defective 
MoSe onto graphene, it is then possible to engineer Van der Waals heterostructures to be 
with new electronic structures determined by DE [20, 21]. 

In this work, we demonstrate that FIB can be used to directly write high-performance 2D 
diodes on a heterostructure consisted of the stacked monolayer MoSe» and graphene. The 
schematic diagram of FIB-writing of the 2D diode on the MoSe;/G heterostructure is 
illustrated in Figure 1. In this prototype, the focused Ga” ion beams are utilized to irradiate the 
MoSe» side of the heterostructure, sputtering Se atoms from the top atomic layer. Se-vacancies 
transfer the irradiated heterostructure from MoSe;-modified graphene to quasi-metallic 
MoSe,/G (x=1.6 or 1.4), according to the fact whether the Dirac point is present. The contact 
of MoSe;,/G-MoSe;G (or MoSe;4/G-MoSe,/G) is similar to the graphene-metal contact, 
generating a p-n junction without the obstacles caused by the Klein tunneling. This work 


offers a candidate route for the direct-writing of electronic circuits on 2D materials via FIB. 


Defect-engineered MoSe,/G heterostructure 


The MoSe,/G heterostructure with MoSe; monolayer on the top is irradiated by a focused 
Ga‘ ion beam with energy of 30 keV and fluence between 10"? and 10 ions/cm’. The 
composition of the heterostructure before and after the FIB-writing is determined through X- 
ray photoelectron spectroscopy (XPS), as shown in Supplementary I. Before writing, the 
heterostructure has stoichiometric ratios of Se:Mo = 2, and C:Mo = 12, suggesting the high 
crystallinity of the as-prepared heterostructure. After irradiation, the integration of the 
photoemission peak shows a reduction in Se-3d, a slight change in Mo-3d, and a negligible 
variation in C-1s atoms. This indicates that numerous Se atoms are removed from the top 
atomic layer, while Mo and C atomic layers remain comparatively unchanged. Figure 2a 
shows the evolution of the Se:Mo stoichiometric ratio (x) as a function of the ion fluence, 
demonstrating stabilized stoichiometries of MoSe, (x = 2, 1.8, 1.6, and 1.4) in agreement with 
previous time-resolved ion sputtering study of Moie: [22-24]. The dangling bonds left by 
sputtered Se atoms are investigated via the comparison of the high-resolution XPS spectra of 
MoSe,/G (Mo-3d) in Fig. 2b and 2c. The Mo* peaks of MoSe»;G locates at 229.2 eV (Mo- 
3d32) and 232.3 eV (Mo-3ds), whilst Mo peaks at ~232.8 eV and ~236.2 eV appear in 
MoSe;4/G, MoSe;./G and MoSe14/G, indicating the generation of dangling bonds of Mo 
atoms due to Se-vacancies [25]. Figure 2c presents the Raman spectra of MoSe,/G 
heterostructures at bands of MoSe, (Aig) and graphene (D, G, and 2D) normalized by the 
maximum intensity of the G band. After the FIB irradiation, a second peak associated with Se- 
vacancies arises adjacent to the Ai, peak [26], whilst there is negligible change in D, G, and 
2D Raman bands, suggesting the well-preserved graphene. 

Figure 2d displays the atomic force microscopy (AFM) and Kelvin probe microscopy 
(KPFM) images of the two inscribed rectangle shapes in the heterostructure corresponding to 
MoSei;;/G and MoSe14/G. In the AFM image, the thickness of the whole heterostructure is 


uniform after the Ga” ion irradiation (Fig. 2e). This means that the amount of Se-vacancies is 
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not large enough to change the morphology, and only the top atomic layer becomes defective. 
A visible boundary between the irradiated and pristine areas is observed in the KPFM image 
(Fig. 2d), and the rectangular imprint becomes even clearer as x decreases. The KPFM image 
represents the work function (associated with surface potential) of the irradiated 
heterostructure ( Wwosewc) altered by Se-vacancies [27], which has a direct relationship with x, 
i.e., 5.2 eV for the MoSe2./G, 5.1 eV for MoSe1,/G, 5.0 eV for MoSei JG, and 4.85 eV for 
MoSe,.4/G, respectively (Fig. 2f and Supplementary II). This circumstance suggests that DE 
does not change the morphology but can easily modify the electrical property of the 
heterostructure. 

To detect the transport characteristics of the defect-engineered heterostructure, we fabricate 
top-gated field-effect transistors (FETs) based on the MoSe,/G. For x equal to 2 and 1.8 (Fig. 
3a and 3c), MoSe,/G based FETs exhibit the typical ambipolar transport behavior with the 
Dirac point shifted from the gate voltage of -2.0 V (x = 2) to 1.5 V (x = 1.8). The observation 
of the Dirac point indicates unique properties of graphene are well preserved in MoSe;/G and 
MoSe;;/G. The shift of the Dirac point means indicates the n-type (x = 2) and p-type (x = 1.8) 
doping of graphene induced by the MoSe, monolayer. The charge carrier mobilities can be 
extracted from the linear regimes on both sides of the current minimum using 
u-(1/C,):(Olsp/0V,), where C, is the back gate capacitance (1x10? F/cm?). Under the same 
measuring conditions, mobilities of MoSe;/G (hole: 17.4 cm?V'!s'; electron: 22 cm’ V's") 
have comparable amounts with the graphene (hole: 41.1 cm?V'!s*; electron: 28 cm? V's’), 
while as the reduction of x, mobilities of MoSe1,/G lower to 0.24 cm? V's” and 0.22 cm? V's". 
The band structures of MoSe,/G are calculated via First-principles as shown in Fig. 3b (x = 2) 
and Fig. 3d (x = 1.8). The band structure of MoSe;/G appears like a simple sum of the isolated 
MoSe; monolayer and graphene. Although Se-vacancies introduce Se-defect-states in the 
bandgap of MoSe; d (Fig. 3d), there is still a Dirac point in the bandgap of MoSe1,/G. This 


demonstrates that these heterostructures display similar properties with graphene near the 


Fermi level, except for a subtle adjustment by the MoSe, monolayer, which has a good 
agreement with FETs results. For MoSe1,/G (Fig. 3e) and MoSe; /G (Fig. S10b) FETs, the V- 
shaped /sp- Va characteristic disappears (replaced with a constant current ungovernable by the 
gate voltage), which can be explained by the calculated band structure in Fig. 3f. The 
characteristic Dirac point vanishes from K, displaced by hyperbolic-like lines. Meanwhile, the 
Fermi level of the heterostructure passes through the band, indicating the quasi-metallic 
behavior of MoSe,/G (x= 1.6 and 1.4). The conductivity (c) of MoSei6/G (MoSe14/G) is 
measured to be 160 Sim (32 S/m), which is far below Au (4.52x10’ S/m). As o is proportional 
to n (Density of electronic states near Fermi surface) [28] (Supplementary X), we can 
conclude n of MoSe; ed (MoSe; AJO) is only a hundred-thousandth (millionth) of Au. 

To summarize, the FIB-writing process generates Se-vacancies with controllable density 
and distribution on the top atomic layer of the MoSe./G. Se-vacancies change the electronic 
properties of this heterostructure, and classify them into MoSe,-modified graphene (x=2 and 
1.8) and quasi-metallic MoSe,/G (x=1.6 and 1.4) according to the presence of the Dirac point. 
Meanwhile, the work function of the heterostructure can be tailored by x from 5.2 eV to 4.85 
eV. The tailored work function allows establishing a charge transfer region (CTR) or potential 
barrier for any charge carrier tunneling at the interface of MoSe,-modified graphene and 
quasi-metal without the obstacles caused by the Klein tunneling, which paves the way to 


produce lateral 2D diodes on the MoSe;/G heterostructure. 


Rectification performance. 


To provide direct evidence for the current suppression, a MoSe»/G heterostructure is half 
irradiated between drain and source electrodes (D and S) to construct a 2D MoSe,/G- 
MoSe,/G rectifying diode. Figure 4a illustrates SEM (scanning electron microscope) and 
KPFM images of the heterostructure and the electrodes (D and S), in which the ultraviolet 
stands for the irradiated region. Different work functions induced by the FIB-writing cause a 


potential barrier between the irradiated and non-irradiated regions, and the carrier tunneling at 
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the barrier is detected through typical |/sp|-Vsp characteristics. Figure 4b shows rectification 
ratios (r = L»ov/L»ov) and |/sp|-Vsp curves for the MoSe,/G-MoSe»/G diodes corresponding to x 
—2.0, 1.8, 1.6 and 1.4 at room temperature. And the corresponding energy band diagrams are 
represented in Fig. 4c. 

The |/sp|-Vsp characteristic of the MoSe;4,/G-MoSe;/G is symmetric. Even there is a 
potential difference (AV = 0.1) at the interface (Fig. 4c) of MoSe;)G and MoSe;;/G. It is 
because both MoSe;,/G and MoSe;;/G have the Dirac point, which enables the low-energy 
quasiparticles in the heterostructure to behave as massless Dirac fermions. According to the 
Klein tunneling effect, massless Dirac fermions have unlimited transport through potential 
barriers. While, the obvious current suppression appears in the MoSe; /G-MoSe»/G (r = 71) 
and is further enhanced with MoSe:4/G-MoSe./G (r = 1.1<10*), where the irradiated 
heterostructure is transferred to the quasi-metal without the Dirac point. The case of the 
MoSe;4/G-MoSe»/G contact is similar to the metal-graphene contact, in which a dipole layer 
is formed by the charge transfer process taking place at the interface (Fig. 4c). The charge 
transfer shifts the Ær of MoSe»/G significantly, which bends the trace of the Dirac point and 
leads to a long screening length due to the small DOS around the Er. Carriers in MoSe;4/G are 
redistributed to screen the work function difference at the interface, and the screen length is 


proportional to ul" 


[28, 29]. Normally, the screen length in metal (like Au) is just a fraction 
of a nanometer due to the large value of n, which can be ignored [29]. But MoSe14/G as a 
quasi-metal has a very low n (six orders of magnitude less than Au). So the screen length of 
MoSe;4/G is hundreds of nanometers, which is non-negligible. The key feature that originates 
from the long screen length in both MoSe14/G and MoSe»/G is the p-n junction that appears 
near the MoSe, 4/G-MoSe,/G contact. 

Further experiments are performed to quantify the stability, filtering, and frequency 
responses of this MoSe; 4/G-MoSe,/G rectifying diode. As displayed in Fig. 4d, performances 


of this diode are stable for the cycle number of 1x10? in air at room temperature. To detect the 
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filtering, MoSe, /G-MoSe;/G rectifying diode is connected to a rectifier circuit with a loading 
resistor of 10 kQ (Fig. 4e), in which the positive/negative input voltage passes/stops with the 
rectification efficiency of 40.6%. Further changing the loading resistance to 40 kQ, the 
rectification efficiency can increase to 75%, indicating the forward resistance of this 
rectifying diode is about 14 kQ. To explore the frequency response, the MoSe,.4/G-MoSe,/G 
filtering diode is connected in parallel with a smoothing capacitor and a loading resistor, and 
simulated with a sinusoidal input. Figure 4f presents the measured voltage output in a rectifier 
configuration, in which the cutoff frequency is around 9.4 MHz. 

Conclusion. 

The FIB-writing has been demonstrated to be a high-efficient and automated way to 
manufacture lateral 2D diodes on graphene-based heterostructures. To demonstrate the 
superiority of this technique, Table 1 compares the performance of the MoSe; /G-MoSe;/G 
diode discussed in this work with previously reported ones. According to the different 
stacking way, 2D diodes can be classified into vertical and lateral heterojunction types. For 
the vertical diodes, devices possess the highest rectification (107—109), with a current density 
of 10°~10° A/cm? at the maximum reverse voltage. These are considered as the most practical 
potential 2D diodes [15, 16]. However, for the 2D lateral diodes realized through chemical or 
electrostatic doping, the rectification is ~10°, and the current density is ^10? A/cm’, thus 
requiring significant improvements. In this work, the MoSe;4/G-MoSe/G diode has a 
rectification of ~10* and a current density of ~6.7 x 10^ A/cm’, representing comparable 


performance to 2D vertical diodes and evident advantages over other lateral diodes. 


Experimental Methods 

Sample preparation 

The graphene and MoSe; monolayers were grown via chemical vapor deposition (CVD). The 
graphene and MoSe, were then transferred to the Si substrate one after the other through the 
wet-chemistry transfer process. A layer of PMMA was first spin-coated onto the substrate 
with MoSe; or graphene and then baked at 100°C for 5 min. Next, it was immersed into a 2- 
mol/L KOH solution for 2 hours to separate the PMMA/graphene (MoSe;) from the substrate, 
and it was swilled fully with DI water. It was then collected on the target substrate and baked 
at 80°C for 2 hours. Finally, the film was immersed in acetone to remove the PMMA. The 
drain and source electrodes were fabricated via electron-beam lithography and evaporation 
(Au/SiO;/Si, Au: 50 nm, S1O;: 50 nm). 

Ga ion beam writing 

The writing of focused Ga' ions was carried out using the FIB technique (model FEI START 
400S). The acceleration voltage was 30 kV, while the beam current was 9.7 pA, 93 pA, and 
0.79 nA. 

Characterization 

The fabricated devices were characterized using an optical microscope (Leica DM2700), and 
SEM (model FEI START 400S with accelerating voltage of 10kV), an AFM, and a KPFM 
(Bruker, Dimension Icon). The electrical and optoelectronic properties of the devices were 
measured via a Keithley 4200A-SCS at room temperature and room pressure. MoSe,/G 
heterostructures used for Raman spectra and XPS spectra measurements were produced by 
CVD and transferred on copper nets. The power of the laser used in Raman spectra 
measurements was 1mW, and the wavelength was 473 nm. 

Top-gated Field-Effect Transistors 

A small droplet of saturated trifluorochloroethylene is dropped onto the top of FETs. The 


probe is inserted into the droplet above MoSe,/G. 
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DFT calculations 

First-principles calculations were performed using the density functional theory (DFT), 
adopting the plane-wave pseudopotential approach as implemented in the Vienna Ab initio 
Simulation Package (VASP) [32]. The electron-electron interaction in the calculations was 
described via the Perdew-Burke-Ernzerhof (PBE) functional under the generalized gradient 
approximation (GGA) [33]. The Van der Waals weak interaction between graphene and 
MoSe» was corrected via the Grimme method (DFT-D2). For the electron wave function, a 
plane wave with a cutoff energy of 500 eV was used. A vacuum region of 30 A was applied in 
the z-direction to avoid interactions between adjacent images. The convergence threshold 
parameters for the optimization and the subsequent calculations were set to be 10° eV 
(energy) and 10? eV/A (force). K-point meshes of 5x5x1 were used in sampling the Brillouin 
zone for the 3x3-MoSe,/4x4-graphene supercell of the heterostructure. All the atomic 
positions and lattice vectors were fully optimized using the conjugate gradient (CG) scheme 


without any restrictions. 
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Figure 1. Schematic diagram of the FIB-writing process. 
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Figure 2. Defect-engineered MoSe,/G heterostructure. (a) Evolution of the stoichiometric 


ratio Se:Mo as a function of the ion fluence of the Ga" ion beam. (b) XPS spectra of MoSe,/G 


at Mo-3d. (c) Raman spectra of MoSe,/G. (Aig peak for MoSe,; D, G, and 2D peaks for 


graphene) (d) AFM and KPFM images of the heterostructures written in the rectangular 


shapes corresponding to MoSe;;/G and MoSe,J/G. 


(e) Thickness/work-function 


corresponding to the pink/purple dashed lines in (d). (f) The measured work function of 


MoSe,/G. 
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Figure 3. Calculated electronic structure of the MoSe,/G heterostructure. Gate dependence 
(Va) of the drain current Isp-Vsp of MoSe2/G (a), MoSe1.s/G (c), and MoSe; JG (e). Calculated 


electronic structures of MoSe;/G (b), MoSe; sl? (d), and MoSe:s6/G (f). The Fermi level is set to 0. 
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Figure 4. Rectifying diodes based on the MoSe,/G heterostructure. (a) SEM (top) and KPFM 
(bottom) images of the MoSe.,/G rectifying diode. (b) Evolution of the rectification ratio as a 
function of x. The inset image is |/sp|-Vsp characteristics of the rectifying diodes constructed with 
MoSe;/G, MoSe18/G, MoSe:«/G, and MoSe14/G, at room temperature. (c) Corresponding energy band 
diagrams of MoSe;3/G-MoSe./G and MoSe;4/G-MoSe./G. (d) The variation ratio of the 
current and rectification (at Vsp = +10 V) of the rectifying diode along with the cycle number. 
(e) Time-resolved input and output voltages of the rectifier circuit. The loading resistor in this 
circuit is 10 kQ. (f) Frequency-resolved output voltage in the rectifier configuration. All 


measurements were taken in air at room temperature. 
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Table 1. Comparison of 2D diodes in literature. 


Heterojunctio 


Laver st : = 
n type ayer structure 7 


1.0x10° at 2 V 21A/cm? at 2 V 
G/MoS; 1.5x10? at 4 V? 3.5x10? A/cm’? at 4 V? [35] 
2.0x10?at 2 V 1.0Ax10°/cm’ at 2 V 


Vertical diodes (ema. 2.0x105 at 5 V° 2.0x10? A/cm2at 5 V° SS 
1.0x10? at 2 V 20 A/cm? at 2 V 


WSejMoS; 1.0x10? at 2 V? 50 A/cm^ at 2 V? [37] 


BP'/MoS; 1.0x10? at 2 V? 1 A/em? at 2 V? 


WSe;MoS; 1.0x10? at 2 V* 33 A/cm? at 2 V° 


1.0x10?at 5 V? 0.7 A/cm? at 5 V? 
WSe/ WSe; 10at2 V 0.02A/cm?at2 V [að] 


Lateral diodes 1.0x10? at 5 V? 5 A/cm’ at 5 V? 


WSejMoSe; -5at2V 0.15A/cm?^at 2 V 


1.1x10^ at 20V ° 2.4x10* A/cm’ at 20V? 


Moses 1.0%10? at 2V 6.510" A/cm? at 2V 


This work 


* The maximum reverse voltage 
^ Black phosphorus 
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